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From a viewpoint of oblate-prolate symmetry and its breaking, we adopt the quadrupole 
collective Hamiltonian to study dynamics of triaxial deformation in shape coexistence phe- 
nomena. It accommodates the axially symmetric rotor model, the 7-unstable model, the 
^ 1 rigid triaxial rotor model and an ideal situation for the oblate-prolate shape coexistence as 

particular cases. Numerical solutions of this model yield a number of interesting suggestions. 
(1) The relative energy of the excited + state can be a signature of the potential shape along 
the 7 direction. (2) Specific E2 transition probabilities are sensitive to the breaking of the 
CO ' oblate-prolate symmetry. (3) Nuclear rotation may induce the localization of collective wave 

functions in the (/3, 7) deformation space. 

o : 

\ §1. Introduction 

! In recent years, experimental data suggesting coexistence of the ground-state 

rotational band with the oblate shape and the excited band with the prolate shape 
have been obtained in proton-rich unstable nuclei. 1 ) -3 ) Together with a variety of 
shape coexistence phenomena observed in various regions of nuclear chart, 4 )~ 7 ) some 
of which involve the spherical shape also, these discoveries stimulate development 
^ ■ of nuclear structure theory capable of describing this new class of phenomena. 8 )~ 13 ) 

Recently, Hinohara et al. 14 ^' 15 ^ carried out detailed microscopic calculations for 
the oblate-prolate shape mixing by means of the adiabatic self-consistent collective 
coordinate (ASCC) method. 8 ) They suggest that the excitation spectrum of 68 Se 
may be regarded as a case corresponding to an intermediate situation between the 
well-developed oblate-prolate shape coexistence limit, where the shapes of the two 
coexisting rotational bands are well localized in the (/?, 7) deformation space, and 
the 7-unstable limit, where a large-amplitude shape fluctuation takes place in the 
7 degree of freedom. Here, (5 and 7 are well-known dynamical variables denoting 
the magnitude of the quadrupole deformation and the degree of axial symmetry 
breaking, respectively. These calculations indicate the importance of the coupled 
motion of the large-amplitude shape fluctuation in the 7 degree of freedom and the 
three-dimensional rotation associated with the triaxial shape. 
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In order to discuss the oblate-prolate shape coexistence phenomena in a wider 
perspective including their relations to other classes of low-energy spectra of nuclei, 
we introduce in this paper a simple phenomenological model capable of describ- 
ing the coupled motion of the large-amplitude 7-vibrational motions and the three- 
dimensional rotational motions. We call it (1+3)D model in order to explicitly 
indicate the numbers of vibrational (7) and rotational (three Euler angles) degrees 
of freedom. This model is able to describe several interesting limits in a unified 
perspective. It includes the axially symmetric rotor model, the 7-unstable model, 16 ) 
the triaxial rigid rotor model 17 ) and an ideal situation of the oblate-prolate shape 
coexistence. It also enables us to describe intermediate situations between these 
different limits by varying a few parameters. This investigation will provide a new 
insight concerning connections between microscopic descriptions of oblate-prolate 
shape mixing and various macroscopic pictures on low-energy spectra in terms of 
phenomenological models. It is intended to be complementary to the microscopic 
approach we are developing on the basis of the ASCC method. 

The (1+3)D model is introduced on the basis of the well-known five-dimensional 
(5D) quadrupole collective Hamiltonian 18 ) by fixing the axial deformation parameter 
/3. Simple functional forms are assumed for the collective potential and the collective 
mass (inertial function) with respect to the triaxial deformation 7. We analyze prop- 
erties of excitation spectra, quadrupole moments and transition probabilities from a 
viewpoint of oblate-prolate symmetry and its breaking varying a few parameters char- 
acterizing the collective potential and the collective mass. Specifically, we investigate 
the sensitivity of these properties to 1) the barrier-height between the oblate and 
prolate local minima in the collective potential, 2) the asymmetry parameter that 
controls the degree of oblate-prolate symmetry breaking in the collective potential, 
and 3) the mass-asymmetry parameter introducing the oblate-prolate asymmetry 
in the vibrational and rotational collective mass functions. Dynamical mechanism 
determining localization and derealization of collective wave functions in the defor- 
mation space is investigated. We find a number of interesting features which have 
received little attention until now: 1) the unique behavior of the excited + state 
as a function of the barrier-height parameter, 2) specific E2 transition probabilities 
sensitive to the degree of oblate-prolate symmetry breaking, 3) the rotation-assisted 
localization of collective wave functions in the deformation space. We also examine 
the validity of the (1+3)D model by taking into account the j3 degree of freedom, 
i.e., by solving the collective Schrddinger equation for a 5D quadrupole collective 
Hamiltonian that simulates the situation under consideration. 

In the present paper, therefore, we intend to clarify the role of the (3-j depen- 
dence of the collective mass. Since the original papers by Bohr and Mottelson, 19 )' 20 ) 
the collective Schrddinger equation with the 5D quadrupole collective Hamiltonian 18 ) 
has been widely used as the basic framework to investigate low-frequency quadrupole 
modes of excitation in nuclei. The collective potential and the collective masses ap- 
pearing in the Hamiltonian have been introduced either phenomenologically 21 ) -24 ) 
or through microscopic calculations. 25 ) _29 )> 31 )- 35 ) in recent years, powerful algebraic 
methods of solving the collective Schrddinger equation have been develop ed 36 )' 37 ) 
and analytical solutions have been found 38 ) -42 ) for some special forms of the collec- 
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tive potential (see the recent review 43 ) for an extensive list of references). However, 
all collective masses are assumed to be equal and a constant in these papers. 38 )~ 42 ) 
It should be emphasized that this approximation is justified only for harmonic vibra- 
tions about the spherical shape. The collective masses express inertia of vibrational 
and rotational motions, so that they play crucial roles in determining the collec- 
tive dynamics. In general, they are coordinate-dependent, i.e., functions of (3 and 
7. In fact, various microscopic calculations for the collective masses indicate their 
significant variations as functions of (3 and ^. 28 )> 29 )> 31 )~ 35 )> 44 )~ 48 ) i n phenomenolog- 
ical analysis of experimental data, for instance, Jolos and Brentano 49 ) have shown 
that it is necessary to use different collective masses for the rotational and (3- and 7- 
vibrational modes in order to describe interband E2 transitions in prolately deformed 
nuclei. 

This paper is organized as follows. After recapitulating the 5D quadrupole 
collective Hamiltonian and the collective Schrodinger equation in §2, we introduce 
in §3 the (1+3)D model of triaxial deformation dynamics. In §4, using the (1+3)D 
model, we first discuss properties of excitation spectra for the collective potentials 
possessing oblate-prolate symmetry. We then investigate how they change when this 
symmetry is broken in the collective potential and/or in the collective mass. In §5 
we examine the validity of the (1+3)D model by introducing (3-j coupling effects 
on the basis of the 5D quadrupole collective Hamiltonian. Concluding remarks are 
given in §6. 

§2. Five-dimensional quadrupole collective Hamiltonian 

We start with the 5D quadrupole collective Hamiltonian involving five collective 
coordinates, i.e., two deformation variables (/?, 7) and three Euler angles: 

H = T vib + T rot + y(/3, 7 ), (2-1) 

T vib = ^Dpp(/3, 7)^ 2 + Dpy(@, 7)^7 + ^Z> 77 (/3, 7)-y 2 , (2-2) 

3 1 

T rot = ^-J fc (/3,7)^. (2-3) 
fc=i 

Here, the collective potential V(/3, 7) is a function of two deformation coordinates, 
(3 and 7, which represent the magnitudes of quadrupole deformation and triaxiality, 
respectively. It must be a scalar under rotation, so that it can be written as a 
function of 1 and /3 3 cos 37. 18 - 1 As is well known, one can restrict the range of 7 to 
be 0° < 7 < 60° by virtue of the transformation properties between different choices 
of the principal axes. The first term T v ;b in Eq. (2-1) represents the kinetic energies 
of shape vibrations; it is a function of (3 and 7 as well as their time derivatives (3 and 
7. The second term T rot represents the rotational energy written in terms of three 
angular velocities ojk which are related to the time derivatives of the Euler angles. 
The three moments of inertia can be written as 



J fc (/?,7) =4/3 2 J D fc (/3,7)sin 2 7fc , 



(2-4) 
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with respect to the principal axes (k = 1 — 3), where 7fc = 7 — (2irk)/3. In this 
paper, we adopt Bohr-Mottelson's notation 18 ) for the six collective mass functions, 
Dpp, -D/3 7 , Dy-y, Di, Z?2 and -D3. They must fulfill the following conditions, 



L> 1 (/3, 7 = o )=L» 2 (A7 = 0°), 

J D 3 (/3,7 = 0°) = J D 77 (/3,7 = )/?" 2 , 
Di(A7 = 60°) = J D 3 (/3,7 = 60 o ), 

D 2 (/3, 7 = 60°) = D 77 (/3,7 = 60°)/T 2 , 



(2-5) 
(2-6) 
(2-7) 
(2-8) 



in the prolate (7 = 0°) and the oblate (7 = 60°) axially symmetric limits. 27 ) 

The classical collective Hamiltonian (2-1) is quantized according to the Pauli 
prescription. Then, the explicit expressions for the vibrational and rotational kinetic 
energies are given by 26 ) 



r v ib 



-h 2 f 1 



2VWR \P 4 



dpi? 





Dp^dry 



p 1 sin 37 



-d^ ( J ^ sin 3ryDp y dp j + <9 7 ( J ^ sin S'yDppd^ 



and 



L rot 



f 2 



fc -2j7- fe (A7)' 
respectively, where and R are the abbreviations of 



W(p,i) = (3~ 2 [Dppi^D^j) - D%(/3,j)] 



(2-9) 



(2-10) 



(2-11) 
(2-12) 



and Ifc are the angular momentum operators with respect to the principal-axis frame 
associated with a rotating deformed nucleus (the body-fixed PA frame). 
The collective Schrodinger equation is written as 



[f vih + r rot + v((3,i)]v IMa (P,i, n) = E Ita & IM *((3,7, fi), 



(2-13) 



where the collective wave function ^iMa{Pil^) is specified by the total angular 
momentum /, its projection M onto the z-axis in the laboratory frame and a dis- 
tinguishing eigenstates possessing the same values of / and M. In Eq. (2-13), Q 
denotes a set of the three Euler angles, which are here dynamical variables describ- 
ing the directions of the body-fixed PA frame with respect to the laboratory frame. 
Using the rotational wave functions T> MK (Q), the orthonormalized collective wave 
functions in the laboratory frame can be written as 



K 



(2-14) 
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where 

(f2\IMK) = + H^m-kM) . (2-15) 

The functions ^iKaiP,^) are vibrational wave functions orthonormalized by 

! dT'J2$*IK a (P,l)$IKa'(P,l) = <W (2-16) 
J K 

with the intrinsic volume element 



dr' = 2^^/W{j3,i)R({5,i) sin 3 7 d/3d 7 . (2-17) 

In Eqs. (2-14) and (2-16), the sum is taken over even values of K from to / for even 
/ (from 2 to / — 1 for odd /). Detailed discussions on the symmetries and boundary 
conditions in the vibrational wave functions <&iKa((3,7) are given, e.g., in Ref. 27). 

§3. Reduction to the (l+3)-dimensional collective Hamiltonian 



For the reason mentioned in §1, we are particularly interested in triaxial defor- 
mation dynamics. In order to concentrate on the 7 degree of freedom, we introduce 
a simple (1+3)D model involving only one vibrational coordinate 7 and three ro- 
tational coordinates. This is done by freezing the (5 degree of freedom in the 5D 
quadrupole collective Hamiltonian (2.1) as explained below. 

The collective potential of this model takes a very simple form: 

y( 7 ) = V sin 2 37 + Vi cos 37. (3-1) 

This form is readily obtained by retaining up to the second order with respect to 
/3 3 cos 37 and fixing f3 at a constant value in an expansion of V(f3,"y) in powers 
of j3 2 and /3 3 cos 37. 50 ) When V\ = 0, the collective potential is symmetric about 
7 = 30° with respect to the transformation 7 — > 60° — 7. For brevity, let us call 
this symmetry and transformation OP (oblate-prolate) symmetry and OP inversion, 
respectively. For positive Vq, two degenerate minima appear at the oblate (7 = 60°) 
and the prolate (7 = 0°) shapes, and they are separated by a barrier that takes the 
maximum at 7 = 30°. We therefore call Vo barrier-height parameter. On the other 
hand, the maximally triaxial shape at 7 = 30° becomes the minimum for negative Vo, 
and it becomes deeper as \Vo\ increases. When V\ 7^ 0, the OP symmetry is broken, 
and the oblate (prolate) shape becomes the minimum for a combination of positive 
Vq and positive (negative) V\. We call V\ asymmetry parameter after its controlling 
the magnitude of the OP symmetry breaking. Thus, by varying the two parameters 
Vq and V\, we can see how the excitation spectrum depends on the barrier height 
and symmetry breaking between the two local minima in the collective potential. 

We present in Fig. 1 some examples of the collective potential V( 7 ) that seem to 
be relevant to a variety of oblate-prolate shape coexistence phenomena. In this figure, 
we can clearly see how the asymmetry between the oblate and prolate minima grows 
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Fig. 1. The collective potentials V(fy) denned by Eq. (2.4) are plotted for the Vq > case with 
solid, dash-dotted and dashed lines as functions of 7 for V1/V0 = 0.0,0.5 and 1.0, respectively. 



as a function of V\ and how the barrier height measured from the second minimum 
sensitively depends on the ratio of V\ to Vq. 

The vibrational kinetic energy term reduces, in the (1+3)D model Hamiltonian, 
to the following form: 

f vib = ~ h2 =^—d, f J ^ ' 7 \ sin3 7 c0 . (3-2) 

We parametrize the collective mass functions as 

D 77 ( 7 ) = D 77 (/3 ,7) = Dffi(l + ecos3 7 ), (3-3) 
D k fr) = D k (J3o,-y) = J D(l + ecos 7fc ). (3-4) 

These are the most simple forms involving only one parameter e that controls the 
degree of OP symmetry breaking in these four mass functions under the requirement 
that they should fulfill the conditions (2-5)— (2-8). We call e mass- asymmetry pa- 
rameter. These functional forms are obtainable also by taking the lowest-order term 
that brings about the OP symmetry breaking in the expressions of the collective 
mass functions microscopically derived by Yamada 48 ^ using the SCC method. 51 ) In 
this paper, we use (3q = 0.1 and D = 50 MeV" 1 , which roughly simulate the values 
obtained in the microscopic ASCC calculation. 14 )> 15 ) 

We solve the collective Schrodinger equation (2-13) replacing V(/?, 7 ), T vi b and 
i7fc(/3, 7) in T rot by Eqs. (3T), (3-2) and i7fc(A)>7)> respectively, with Eqs. (3-3) and 
(3-4). Accordingly, the collective wave function is denoted &iMa{l,ft)- Note that 
the sign change V\ — > —V\ corresponds to the OP inversion. One can then easily 
confirm that the OP inversion is equivalent to the simultaneous sign change of the 
parameters, (e,±Vi) — ► (— e, T^i)) m the (1+3)D model Hamiltonian. Therefore, it 
is enough to study only the case of positive e. 

Figure 2(a) shows, for an example of e = 0.5, how the collective mass functions 
-D 77 and behave as functions of 7. Figure 2(b) indicates the degree of oblate- 
prolate asymmetry in the moments of inertia Jk{l) brought about by the terms 
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Fig. 2. (a) The vibrational and rotational collective mass functions D 77 (7) and £^(7), defined by 
Eqs. (2.2) and (2.3), are plotted as functions of 7 for the mass-asymmetry parameter e = 0.5. 
(b) Three rotational moments of inertia •Jk(f3,'y) with k =1, 2 and 3, denned by Eqs. (1.2) and 
(2.3), are plotted with solid lines as functions of 7 for the mass-asymmetry parameter e = 0.5. 
The magnitude of the quadrupole deformation j3 is fixed at /3 2 = 0.1. For comparison, the 
moments of inertia for the t = case are also plotted with dashed lines. 

involving e. From these figures we can anticipate that, for positive (negative) e and 
under the parameterizations of Eqs. (3-3) and (3-4), the rotational energy favors 
the oblate (prolate) shape while the vibrational energy prefers the prolate (oblate) 
shape. 

§4. Triaxial deformation dynamics 

We have built a new computer code to solve the collective Schrodinger equa- 
tion (2-13) for general 5D quadrupole collective Hamiltonian as well as its reduced 
version for the (1+3)D model. Numerical algorithm similar to that of Kumar and 
Baranger 27 ) is adopted in this code, except that we discretize the (/?, 7) plane us- 
ing meshes in the (5 and 7 directions in place of their triangular mesh. Numerical 
accuracy and convergence were checked by comparing our numerical results with ana- 
lytical solutions in the spherical harmonic vibration limit, in the so-called square- well 
type /3 dependence limit and in the 7- unstable model. 16 ) In the following, we discuss 
how the solutions of the collective Schrodinger equation depend on the barrier-height 
parameter Vo, the asymmetry parameter V\ and the mass- asymmetry parameter e. 
We use the adjectives yrast and yrare for the lowest and the second-lowest states for 
a given angular momentum I, respectively, and distinguish the two by suffices as I\ 
and I2. 

4.1. Excitation spectra in the presence of OP symmetry 

Let us first discuss the situation where V± = and e = 0. In this case, both the 
collective potential and the collective mass functions are symmetric about 7 = 30°, 
so that the (1+3)D model Hamiltonian possesses OP symmetry. Furthermore, the 
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collective mass parameter D and (3q enter the collective Schodinger equation only in 
the form of the overall factor (2D/3q) _1 in the kinetic energy terms. Hence, only the 
ratio of the barrier-height parameter Vq to this factor is important to determine the 
collective dynamics. In the particular case of Vo = 0, i.e., V{pf) = 0, which is well 
known as the Wilets-Jean 7-unstable model, 16 ) the excitation spectra are completely 
scaled by this factor. 

In Fig. 3, we show the excitation spectra for the case of V\ = and e = as 
functions of Vo- Here, the excitation energies are normalized by the excitation energy 
of the first excited 0+ state (denoted Oj) at V = 0, E(Q%) (which is 1.8 MeV for 
Pq = 0.1 and D = 50 MeV -1 adopted in the present calculation). Accordingly, this 
figure is valid for any value of (2D/3g) _1 by virtue of the scaling property mentioned 
above. In the lower panels of this figure, the collective potentials V(j) and the 
ground + state energies E(0J) are illustrated for three typical situations: 1) a 




10 20 30 40 50 60 10 20 30 40 50 60 -20 20 40 60 80 



y y y 

(b) (c) (d) 

Fig. 3. In the upper panel (a), dependence of excitation spectrum on the parameter Vo is displayed. 
Excitation energies as well as Vo values are normalized by the excitation energy E{0^ ) of the 
second (the first excited) + state for V(j) = 0. In the lower panels (b), (c) and (d), the 
potentials V(j) and the ground + state energies E(0{) are illustrated for three different values 
of Vo/E(0%) = -5.0,0.0 and 5.0, respectively. Note that E(0f) = for Vo = 0. Because the 
collective potential V(y) is a periodic function of 60° in 7, only the region 0° < 7 < 60° is 
drawn with a solid line in (d). 
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triaxially deformed case where a deep minimum appears at the triaxial shape with 
7 = 30°, 2) the 7-unstable case, where the collective potential is flat with respect to 
7, and 3) an extreme case of shape coexistence where the oblate and prolate minima 
are exactly degenerate in energy. (Strictly, shape coexistence does not appear in a 
case where the two minima are exactly degenerate as we shall see in Fig. 9.) Note 
that the collective potential V{~j) is a periodic function of 60° in 7. Therefore, it is 
displayed by the solid line only in the region 0° < 7 < 60°. 

In the positive- Vo side of this figure, it is clearly seen that a doublet structure 
emerges when the barrier-height parameter Vq becomes very large. In other words, 
approximately degenerate pairs of eigenstates appear for every angular momentum 
when Vo/i?(0^) ^ 1. This is nothing but the doublet pattern known well in the 
problems of double-well potential. 52 ) In the present case, this doublet structure is 
associated with the OP symmetry. Furthermore, one immediately notices a very 
unique behavior of the Oij" state. When the barrier-height parameter Vo decreases 
(from the limiting situation mentioned above), its energy rises more rapidly than 
those of the yrare 2^~, and states. Thus, a level crossing of the and 
states takes place at Vo/-E(0^") ~ 3. When Vo further decreases and approaches zero, 
the excitation energy of the 0^ state approaches those of the 4^ and states. In 
the 7-unstable limit of Vo = 0, they are degenerate. 

In the negative- Vb side, the excitation energies of the 3^~ and 5^ states drastically 
decrease as Vo decreases, and the excitation spectrum characteristic to the Davydov- 
Filippov rigid triaxial rotor model 17 ) appears when the triaxial minimum becomes 
very deep, i.e., when Vq/E^) "C — 1. 
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4.2. Breaking of the OP symmetry in the collective potential 

Next, let us investigate effects of the OP-symmetry-breaking term 14 cos 37 in 
the collective potential V(j). The effects are manifestly seen for the case of e = 
and Vq = 0. We present in Fig. 4 the excitation spectrum for this case as a function 
of the asymmetry parameter V\. For V\ = 0, the spectrum exhibits the degeneracy 
characteristic to the 7-unstable model: 16 ) e.g., the 1^1 1^2 an< ^ ®t states are 
degenerate. With increasing the magnitude of V\, such degeneracies are lifted and 
the yrare O2 ,2^" ^2 an< ^ ^2 states as well as the odd angular momentum 3^~ and 

states rise in energy. As a consequence, the well-known ground-state rotational 
band spectrum appears for sufficiently large values of |Vi|. For instance, we can see in 
Fig. 4 that, as V\ increases, the ratio of the yrast 2^~ and 4^ energies, E(4f)/E(2f), 
increases from 2.5 at V\ = 0, which is the value peculiar to the 7-unstable model, 
to be 3.3 for sufficiently large V\. This figure beautifully demonstrates the fact that 
the breaking of spherical symmetry is not sufficient for the appearance of regular 
rotational spectra even if the magnitude of quadrupole deformation is considerably 
large: we also need appreciable amount of the OP symmetry breaking. We also note 
that the spectrum does not depend on the sign of V\. This is because, for e = 0, 
both the vibrational and rotational kinetic energy terms in the (1+3)D collective 
Hamiltonian possess the symmetry under the OP inversion. In short, the inversion 
V\ — > —V\ merely interchanges the roles of the oblate and the prolate shapes. 
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Let us proceed to more general situations where e = but both Vq and V\ are 
nonzero. Three typical situations are illustrated in Fig. 5. In this figure, the collec- 
tive potentials and energy spectra are drawn for three different sets of parameters, 
(Vb,Vi) = (— 2, 1), (1, 0.5) and (10, 1) MeV. Panel (a) simulates a situation where 
the minimum of the collective potential occurs at a triaxial shape but it is rather 
shallow, so that the potential pocket accommodates only the ground + state and 
the first excited 2 + state. Panel (b) simulates the situation encountered in the mi- 
croscopic ASCC calculation, 14 ^' 15 ) where two local minima appear both at the oblate 
and prolate shapes but the barrier between them is so low that strong shape mixing 
may take place. Panel (c) illustrates an ideal situation for shape coexistence, where 
the barrier between the oblate and prolate minima is so high that two rotational 
bands associated with them retain their identities. 

For three typical situations illustrated in Fig. 5, we examine the sensitivity of 
excitation spectra to the asymmetry parameter V\. This is done in Fig. 6, (Panels 
(a), (b) and (c) display the excitation spectra as functions of V\ for Vq = —2, 1 and 
10 MeV, respectively.) We see that the dependence on V\ is rather weak in the cases 



> 

CD 



CT> 

CD 
£Z 
LU 

tz 
o 

s 

lo 

LU 



















ft- 1 

■ 3^ 






4; 


- 2| 







3.5 
3 
2.5 



>, 
en 

05 2 

LU 

c 1-5 
o 

§ 1 
'o 





1 


1 








6 i : 




















— - 








5 : 


































— * 






2 i ; 


-61 




— ' 








4;- 


^ — ~~ 


























2;- 


1 1 1 1 



0.5 1 

V 1 [MeV] 

(a) V = -2.0MeV 

3.5 

5 



1.5 



0.2 0.4 0.6 
V1 [MeV] 

(c) Vo = lO.OMeV 



> 

CD 



5*2.5 



CD 2 
LU 

c 1.5 
o 

S 1 

o 

x 0.5 

LU 



i 



,4, 



67" 



^2: 




0.2 0.4 0.6 
V 1 [MeV] 

(b) Vo = l.OMeV 



0.8 



0.8 



Fig. 6. Dependence of excitation spectrum on the asymmetry parameter Vi is displayed for Vo = 
-2.0 (panel a), 1.0 (panel b) and 10.0 MeV (panel c). 
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(a) and (b). In contrast, the effect of the V\ term is extremely strong in the case 
(c): the approximate doublet structure at Vo = is quickly broken as soon as the V\ 
term is switched on, and the excitation energies of the yrare partners (0^~, 2^~, 4^, 6^~) 
remarkably increase as V\ increases. Quantitatively, the energy splittings between 
the yrast and yrare states with the same angular momenta are proportional to 
f° r Vi 0.2 MeV and then increase almost linearly in the region of V\ > 0.2 MeV. 
The quadratic dependence on V\ in the small- V\ region can be understood as the 
second-order perturbation effects in the double-well problem. In this case, what is 
important is not the absolute magnitude of V\ but its ratio to the energy splitting due 
to the quantum tunneling through the potential barrier between the two minima. 
The structure of the collective wave functions is drastically changed by the small 
perturbation. Indeed, for V\ > 0.2 MeV, they are already well localized in one of 
the potential pockets, as we shall discuss below in connection with Fig. 9. Once the 
collective wave functions are well localized, the energy shifts are mainly determined 
by the diagonal matrix elements of the V\ term (with respect to the localized wave 
functions), leading to the almost linear dependence on V\. 

The effects of the V\ term are seen more clearly in the collective wave functions. 
Figures 7, 8 and 9 show the collective wave functions squared YIk \^iKa{l)\ 2 of the 
yrast and yrare states for Vq = —2, 1 and 10 MeV, which respectively correspond to 
the potentials in panels (a), (b) and (c) of Fig. 5. In Figs. 7 and 8, one sees that the 
localization (with respect to the 7 coordinate) of the collective wave functions of the 
yrast states grows as V\ increases, while that of the yrare states are insensitive to 
V\. This can be easily understood by considering the kinetic-energy effects tend to 
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Fig. 7. Collective wave functions squared J2 K |^/ft- Q (7)| 2 of the yrast and the yrare states for the 
collective potentials V(7) with Vi = 0.0, 1.0 and 2.0 MeV. The barrier-height parameter Vo is 
fixed at -2.0 MeV. 
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dominate over the potential-energy effects with increase in the excitation energies. 
It is interesting to note that in Fig. 8, although the yrast 4^~ state is situated above 
the potential barrier, its wave function is well localized around the oblate minimum. 
Concerning the yrare state, although its wave function has the maximum at the 
prolate shape due to the orthogonality to the 4^ state, it is considerably extended 
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over the entire region of 7. For the reason mentioned above, the effects of the V\ 
term on localization properties of the yrare states are much weaker than those of the 
yrast states. 

In Fig. 9, we see drastic effects of the V\ term. At V\ = 0, the wave functions 
squared are symmetric about 7 = 30°. This symmetry is immediately broken after 
the V\ term is switched on: the wave functions of the yrast states rapidly localize 
around the potential minimum even if the energy difference between the two local 
minima is very small. In striking contrast to the situation presented in Fig. 8, the 
wave functions of the yrare states also localize about the second minimum of the 
potential. We confirmed that V\ of 0.2 MeV is sufficient to bring about such strong 
localization. As pointed out above in connection with the double-well problem, this 
value of V\ is small but comparable to the energy splittings. Thus, one can regard 
Fig. 9 as a very good example demonstrating that even small symmetry breaking in 
the collective potential is able to cause a dramatic change in the properties of the 
collective wave function, provided that the barrier is sufficiently high. 

In the remaining part of this subsection, we discuss how the V\ term affects 
properties of electric quadrupole (E2) transitions and moments. In Fig. 10, the 
B{E2) ratios that vanish in the V\ = limit are plotted as functions of V\ for the 
Vo = —2.0 MeV case corresponding to Figs. 5(a), 6(a) and 7. Vanishing of these 
ratios is well known as one of the signatures of the triaxial shape with 7 = 30° 
in the rigid triaxial rotor model. 17 -* It should be emphasized, however, that this 
is in fact a consequence of OP symmetry: these E2 transitions vanish due to the 
exact cancellation between the contribution from the prolate side (0° < 7 < 30°) 
and that from the oblate side (30° < 7 < 60°). Therefore, the localization around 
7 = 30° is not a necessary condition. In fact, these ratios vanish also in the 7-unstable 
model. 16 ) As anticipated, we see in Fig. 10 that these ratios increase as the collective 
potential becomes more asymmetric with respect to the oblate and prolate shapes. 
In particular, the significant rises of the ratios, S(£'2;4^" — > 3f)/ B(E2; 4^ — > 4^) 
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Fig. 10. Dependence on the asymmetry parameter Vi of selected E2 branching ratios that vanish 
in the oblate-prolate symmetric limit V\ = 0. The barrier-height parameter Vo is fixed at —2.0 
MeV. 
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Fig. 11. Dependence on the asymmetry parameter Vi of the branching ratios B(E2; 7 yrar e — > 7 yra st)/ 
B(E2; Jyraro — > (7 — 2) yrare ). Here, 7 yr ast and 7 yr are denote the yrast and the yrare states having 
the same values of angular momenta 7, while (7— 2) yrarc the yrare states with angular momentum 
7 — 2. The left and right panels display the results of calculation for Vo = 1.0 and 10.0 MeV, 
respectively. 



and B(E2; 6% -> 5f)/B(E2; 6J -> 6J"), are remarkable. This may be interpreted 
as an incipient trend that the sequence of states (2^, 3^", 4^~, 5^ and 6^) forms a 
7-vibrational bandlike structure about the axially symmetric shape (oblate in the 
present case) when V± becomes very large. 

In Fig. 11, the ratios B(E2; I yrarc ->■ I yrast )/ B(E2; J yrare ->■(/- 2) yrarc ) are 
plotted as functions of V\ for the two cases of Vq = 1.0 and 10.0 MeV. The Vq = 1.0 
case corresponds to Figs. 5(b), 6(b) and 8, while the Vq = 10.0 case corresponds to 
Figs. 5(c), 6(c) and 9. Here, the numerator denotes B(E2) values for Al = 0, E2 
transitions from the yrare to the yrast states with the same angular momenta /, 
while the denominator indicates Al = —2, E2 transitions between the yrare states. 
One immediately notices a sharp contrast between the two cases: when the barrier 
between the oblate and prolate local minima is very low (Vo = 1-0 MeV), these values 
are sizable even at V\ = l.OMeV, although they gradually decrease with increase of 
V\. In contrast, when the barrier is very high (Vq = 10.0 MeV), they quickly decrease 
once the OP-symmetry-breaking term is turned on. The reason why the interband 
E2 transitions almost vanish is apparent from Fig. 9; the collective wave functions 
of the yrast and yrare states are well localized around the oblate and prolate shapes, 
respectively. 

One can further confirm the same point by looking at the spectroscopic quadrupole 
moments displayed in Fig. 12 in a way parallel to Fig. 11. They vanish in the pres- 
ence of the OP symmetry due to the exact cancellation between the contributions 
from the oblate side and from the prolate side. In Fig. 12(a) for the low barrier case 
(Vo = 1.0 MeV), the quadrupole moments of both the yrast and yrare states first 
increase after the V\ term is switched on, but those for the yrare states eventually 
saturate. This trend is obvious especially for the and 6^~ states. These results are 
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Fig. 12. Dependence of the spectroscopic quadrupole moments on the asymmetry parameter Vi, 
calculated for Vb = 1.0 MeV (left panel) and Vb = 10.0 MeV (right panel). Their values are 
plotted in units of the intrinsic quadrupole moment, eQo — 3e/\/57rZ7?o/3o- 



easily understandable from the properties of their wave functions displayed in Fig. 8. 
That is, while the localization in the yrast states develops with V\ increasing, the 
wave functions of the yrare states widely extend over the entire region of 7, and the 
effects of the V\ term are rather weak. In contrast, Fig. 12(b) for the high barrier 
case (Vb = 10.0 MeV) demonstrates that both the yrast and yrare states quickly 
acquire quadrupole moments as soon as the V\ term is switched on. This is a direct 
consequence of the wave function localization displayed in Fig. 9. Note again that 
the sign change, V\ <-> —VI, corresponds to the OP inversion. 

When Vb is large and V\ is small, just as in the above case, the yrast and 
the yrare states can be grouped, in a very good approximation, into two rotational 
bands: one associated with the oblate shape and the other with the prolate shape. 
This is an ideal situation for the emergence of an oblate-prolate shape coexistence 
phenomenon. According to the realistic HFB calculations for the collective poten- 
tial, 53 ) however, it seems hard to obtain such a large value of Vq. Therefore, for the 
shape coexistence phenomena we need to take into account dynamical effects going 
beyond the consideration on the collective potential energies. We shall discuss this 
point in the succeeding subsection. 

4.3. Breaking of OP symmetry in the collective mass 

We examine dynamical effects on the localization properties of the collective 
wave functions. As mentioned in §1, we are particularly interested in understanding 
the nature of the shape coexistence phenomena observed in nuclei for which approxi- 
mately degenerate oblate and prolate local minima and a rather low barrier between 
them are suggested from the microscopic potential energy calculations. 14 )' 15 ^' 53 ^ In 
the followings, we therefore concentrate our discussion on the case of the collective 
potential with a low barrier (Vb = 1.0) and weak OP asymmetry (VI = 0.5) repre- 
sented in Fig. 5(b). We shall investigate how the results discussed in the previous 
subsections for the e = case, where the collective mass D 77 (7), -Dfc(7) and the 
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rotational moments of inertia Jk{l) are symmetric functions about 7 = 30°, are 
modified when the mass-asymmetry parameter e becomes nonzero. 

In Fig. 13, the collective wave functions squared |^/Xa(7)| 2 calculated for 
e = 0.5 are displayed. It is clearly seen that, for the yrast states, the localization 
around the oblate shape (7 = 60°) develops with increase in the angular momentum. 
The reason is easily understood: for positive e, the rotational moments of inertia 
perpendicular to the oblate symmetry axis (2nd axis), ^1(7 = 60°) and ^3(7 = 
60°), are larger than those perpendicular to the prolate symmetry axis (3rd axis), 
i7i(7 = 0°) and ^2(7 — 0°), as shown in Fig. 2(b). Therefore, the rotational energy 
for a given angular momentum decreases by increasing the probability of existence 
around the oblate shape. Since the rotational energy dominates over the vibrational 
and potential energies, the localization is enhanced for higher angular momentum 
states. We call this kind of dynamical effect rotation- assisted localization. On the 
other hand, the wave functions of the yrare states exhibit two-peak structure: the 
first peak at the prolate shape (7 = 0°) and the second at the oblate shape(7 = 60°) 
except in the 2^ state. One might naively expect that the yrare states would localize 
about the prolate shape because of the orthogonality requirement to the yrast states. 
However, the second peak is formed around the oblate shape in order to save the 
rotational energy keeping the orthogonality condition. 

Due to the two peak structure of the yrare wave functions mentioned above, the 
ratios of B(E2) values from an yrare state to an yrast state to those between the 
yrare states remain rather large for a wide region of the mass-asymmetry parameter 
e. This is shown in Fig. 14. 

In spite of the two-peak structure of the yrare wave functions, we can find a 
feature of shape coexistence in the spectroscopic quadrupole moments Q, which are 
shown in Fig. 15 as functions of e. Let us first concentrate on the e > part of 
this figure. We see that the Q values of the yrast states are positive, indicating 
their oblatelike character, while the yrare states have negative Q, indicating their 
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Fig. 14. Dependence on the mass-asymmetry parameter e of the ratios B(E2;I yla 



* yrast 



)/B(E2;I y 



(I - 2) yrare ), calculated for Vb = 1.0 and Vi = 0.5 MeV. Here, L 
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and /yrare denote the yrast and the yrare states having the same values of angular momenta I, 
while (I — 2)yrare the yrare states with angular momentum 1—2. 
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Fig. 15. Dependence of the quadrupole moments on the mass-asymmetry parameter e, calculated 
for Vb = 1.0 and Vo = 0.5 MeV. The left and the right panels display the results of calculation 
for the yrast and the yrare states, respectively. 



prolatelike character. This can be regarded as a feature of the shape coexistence. 
It is furthermore seen that the absolute magnitude of Q increases as e increases, 
except the state. As discussed above in connection with Fig. 13, for positive 
e, the oblate shape is favored to lower the rotational energy. Consequently, with e 
increasing, the collective wave functions in the yrast band tend to localize around 
the oblate shape more and more and those of the yrare states localize around the 
prolate shape because of the orthogonality if the angular momentum is not so high. 

It is also noticeable that the absolute magnitude of Q in the yrare band decreases 
with increase in the angular momentum. This is because the cancellation mechanism 
between the contributions from the oblatelike and prolatelike regions of the collective 
wave function works more strongly as the two-peak structure grows. 
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Next, let us discuss on the e < part of Fig. 15. For negative e, the prolate shape 
is favored to lower the rotational energy. On the other hand, the collective potential 
under consideration (Vi = 0.5 MeV) is lower for the oblate shape. Hence, the 
rotational energy and the potential energy compete to localize the collective wave 
function into the opposite directions. It is seen in Fig. 15 that the spectroscopic 
quadrupole moments of the yrast states decrease with e decreasing and that this 
trend is stronger for higher angular momentum states. As a consequence, the Q 
values of the 6^ and 8f states become negative for large negative e, which implies 
that the rotational effect dominates there. 

Finally, we show in Fig. 16 how the localization properties of the collective wave 
functions depend on the mass-asymmetry parameter e, taking the 2 + and 8 + states 
as representatives of low and high angular momentum states. It is clearly seen that, 
for the yrast states, while the the localization of the 2f state around the oblate 





Fig. 16. Dependence of the collective wave functions on the mass-asymmetry parameter e, calcu- 
lated for Vo = 1.0 and Vi = 0.5 MeV. The upper-left and the upper-right panels display the 
results of calculation for the yrast and the yrare states with I n = 2 + , while the lower- left and 
the lower-right panels show those for the yrast and the yrare states with I n = 8 + . 
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shape is rather insensitive to e, that of 8^~ state remarkably develops with increasing 
e. For the yrare states, the localization of the 2Jj~ state around the prolate shape 
grows with increasing e, while the state retains the two-peak structure discussed 
above for every value of e. The different effects of the mass-asymmetry parameter e 
on the 2 + and 8 + states are comprehensible from the consideration of the relative 
importance of the rotational and vibrational energies shown in Table I. We see that 
the rotational energies dominate in the 8 + states, while the vibrational energies are 
comparable in magnitude to the rotational energies in the 2 + states. Thus, the 
effect of the mass-asymmetry parameter e on the localization properties of the 8 + 
states can be easily accounted by the rotational energy. On the other hand, in the 
situation characterized by the parameters Vq = 1.0 and V\ = 0.5 MeV, the properties 
of the 2 + states are determined by a delicate competition between the rotational and 
vibrational kinetic energies as well as the potential energy. The growth of the prolate 
peak with increasing e in the 2ij~ state turns out due mainly to the increase of the 
vibrational mass -D 77 at the prolate shape. 

Table I. Expectation values of the vibrational and rotational energies in units of MeV for the 
2^, 2^~, 8^and 8 J states. The results for e = 0.0 and 0.5 are shown. 
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Summarizing this subsection, we have found that the OP symmetry breaking 
in the collective mass plays an important role in developing the localization of the 
collective wave functions of the yrast states. On the other hand, the asymmetry of 
the collective mass tends to enhance the two peak structure of the yrare states. 

§5. Role of /3-7 couplings 

In this section, we examine whether or not the results obtained in the previous 
section using the (1+3)D model remain valid when we take into account the (3 degree 
of freedom. Because this is a vast subject, we here concentrate on the situation in 
which we are most interested: namely, the case with Vq = 1.0, V\ = 0.5 MeV and 
e = 0.5. 

5.1. A simple (2+ 3) -dimensional model 

We come back to the collective Schrodinger equation (2-13) for the general 5D 
quadrupole collective Hamiltonian and set up the collective potential in the following 
form: 

V{p,i) = ic(/? 2 -/? 2 ) 2 -z;o/3 6 cos 2 37 + z;i/3 3 cos37 + C 6 /? 6 , (5-1) 

where Vq = Vo/(3q and v\ = Vi//3q. The first term ensures that the collective wave 
functions localize around /3 ~ /3q . The second and third terms are reduced to the col- 
lective potential V(j) in the (1+3)D model when the collective coordinate (3 is frozen 
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Fig. 17. Map of the two-dimensional collective potential V(j3, 7) defined by Eq. (5-1) for C = 
800.0, Vb = 1.0, Vi = 0.5, C 6 = 1000.0 MeV and /5g = 0.1. 



at /? = /?o- The fourth term guarantees that the potential satisfies the boundary con- 
dition, V((3, 7) — > 00 as f3 — ► 00. Obviously, the collective potential (5-1) fulfills the 
requirement that it should be a function of /3 2 and /3 3 cos 37. We note that various 
parameterizations of the collective potential similar to Eq.(5T) have been used by 
many authors. 27 ^' 37 )> 40 ) It is certainly interesting and possible to derive the coeffi- 
cients C, vq,vi, Cq and (5q using microscopic theories of nuclear collective motion. In 
this paper, however, we simply treat these coefficients as phenomenological param- 
eters and determine these values so that the resulting collective potential V(f3, 7) 
qualitatively simulates that obtained using the microscopic HFB calculation 53 ) for 
68 Se . They are C = 800, V = 1.0, Vi = 0.5, C 6 = 1000 MeV and ffi = 0.1. In 
Fig. 17, the collective potential with these coefficients is drawn in the (/?, 7) plane. 
One may immediately notice the following characteristic features of this collective 
potential. 1) There are two local minima, one at the oblate shape and the other at 
the prolate shape. They are approximately degenerate in energy but the oblate min- 
imum is slightly lower. 2) There is a valley along the (3 = (3q line connecting the two 
local minima. 3) The spherical shape is a local maximum, which is approximately 4 
MeV higher than the oblate minimum. 

We set up the collective mass functions appearing in the collective kinetic energy 
terms as follows: 



Dp P {p,j) =Z?(l-e / /3cos3 7 ), (5-2) 

D 77 (/3, 7) = D(3 2 (l + e'(3 cos 37), (5-3) 

Dfofai) sin 3 7 , (5-4) 

Dfc(/3,7) =D(l + e7?cos 7fc ), (5-5) 

where e' = e/@Q. These expressions are adopted to take into account the lowest-order 
(/?, 7) dependence of the collective mass functions derived by Yamada 48 ) by means 
of the SCC method. 51 ^ We use the same values for D and e as in Fig. 13, namely 
L>=50MeV" 1 and e = 0.5. 
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5.2. Comparison of the (1+3)D and the (2+3)D model calculations 

We have solved the collective Schrodinger equation (2-13) using the collective po- 
tential (5-1) and the collective masses (5-2)-(5-5). Below, the results of the numerical 
calculation are presented and compared with those of the (1+3)D model. 

Figures 18 and 19 display on the (/?, 7) plane the two-dimensional collective wave 
functions squared, J2k \@iKa{P> 7)| 2 > an d the (3 -weighted ones, fl 4 ^2 K l&iKaiP, 7)| 2 5 
respectively, of the yrast and yrare states with even angular momenta 7 = — 8. 
The P factor carries the major (5 dependence of the intrinsic volume element dr' 
given by Eq. (2-17). We see in Fig. 18 that, except the + states, the yrast (yrare) 
wave functions are well localized around the oblate (prolate) shape. While the local- 
ization of the yrast wave functions grows as the angular momentum increases, the 
yrare wave functions gradually develop the second peaks around the oblate shape. 
These behaviors are qualitatively the same as those we have seen for the collective 
wave functions in the (1+3)D model in Fig. 13. For the + states, the wave func- 
tions squared in Fig. 18 appear to be spread over a rather wide region around the 
spherical shape, but the /3 4 -weighted ones in Fig. 19 take the maxima (as functions 
of (3) near the constant-/? line with (5 = Pq. Thus, we can see in Fig. 19 a very good 
correspondence with the (1+3)D wave functions including the + states also. That 
is, the behaviors along the constant-/? line of the collective wave functions in the 
(2+3)D model exhibit qualitatively the same features as those of the (1+3)D model. 
(A minor difference is seen only in the relative heights of the oblate and prolate peaks 
of the 0^ wave function.) We note that the localization properties of the collective 
wave function are seen better in the /3 4 -weighted wave functions squared than those 
multiplied by the total intrinsic volume element, dr' YIk \^iKa{P, 7)| 2 > which always 
vanish at the oblate and prolate shapes due to the sin 37 factor contained in dr' . 

Finally we compare in Fig. 20 the excitation spectrum obtained in the (2+3)D 
model with that in the (1+3)D model. It is evident that they agree very well. Aside 
from the quantitative difference that the excitation energies are a little higher and 
the B(E2) values increase more rapidly with increase in the angular momentum in 
the (2+3)D model than in the (1+3)D model, the essential features of the excitation 
spectrum are the same. In particular, we find a beautiful agreement in the level 
sequences: the energy ordering of these eigenstates are exactly the same between 
the two calculations. This agreement implies that the /3-7 coupling plays only a sec- 
ondary role here, and the major feature of the excitation spectrum is determined by 
triaxial deformation dynamics. In this dynamics, the 7 dependence of the collective 
mass functions plays an important role as well as that of the collective potential. 

The excitation spectra of Fig. 20 are quite different from any of the patterns 
known well in axially symmetric deformed nuclei, in the rigid triaxial rotor model 
and in the 7-unstable model. It also deviates considerably from the spectrum ex- 
pected in an ideal situation of the oblate-prolate shape coexistence where two rota- 
tional bands keep their identities without strong mixing between them. Among a 
number of interesting features, we first notice the unique character of the 0^~ state. 
It is significantly shifted up in energy from its position expected when the yrare 
0<j~, 2~2 , 4^2, $2 and 8^" states form a regular rotational band. As we have discussed 
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Fig. 18. Collective wave functions squared ^2 K \<&iKa (/?, 7)| 2 calculated for the mass-asymmetry 
parameter e — 0.5 and the two-dimensional collective potential V(/3, 7) with C = 800.0, Vq = 
1.0, Vi = 0.5, C 6 = 1000.0 MeV and ffi = 0.1. The upper and the lower panels display the 
results for the yrast and the yrare states, respectively. 
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Fig. 19. The same as Fig. 18 but for weighted collective wave functions squared 
/3 4 Ea-I^a-«(/3, 7 )| 2 . 



in connection with Fig. 3, the position of the 0^" state relative to the 2jj~ state serves 
as a sensitive measure indicating where the system locates between the 7-unstable 
situation (Vb = V\ = 0) and the ideal oblate-prolate shape coexistence (large Vq and 
small Vi). Thus, the results of our calculation suggest that experimental data for the 
excitation energy of the state provide very valuable information on the barrier 
height between the oblate and prolate local minima. In this connection, we also note 
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Fig. 20. Comparison of the excitation spectra and the _E2-transition properties in the (1+3)D model 
(left side) and those in the (2+3)D model (right side), calculated for the potential parameters 
Vo — 1.0, Vi = 0.5 MeV and the mass-asymmetry parameter e = 0.5. Additional parameters 
for the two-dimensional potential V(0, 7) in the (2+3)D model are C = 800.0 and Ce = 1000.0 
MeV. The B(E2) values are written on the transition arrows normalizing B{E2;2\ — > 0i) as 
100. Weak E2 transitions whose B(E2) values smaller than 10 are not shown. 



that the 3J (5J~) state is situated slightly higher in energy than the 6J and 4^~ (8J 
and 62") states. These are other indicators suggesting that the system is located in 
an intermediate situation between the two limits mentioned above. 

In Fig. 20, we furthermore notice interesting properties of the E2 transitions 
from the yrare to the yrast states: for instance, the E2 transitions with AI = —2, 
or B(E2; I yrarc -»• (J - 2) yrast ), and those with AI = +2, or B(E2; !-►(/ + 2) yrast ), 
are much smaller than other yrare-to-yrast E2 transitions. These transitions are 
forbidden in the 7-unstable model 16 ) because they transfer the boson seniority v 
by Av = 2 and Av = 0, respectively. Thus, some features characteristic to the 
7-unstable situation persist here. On the other hand, although B(E2;3\ — > ^2) 
and B{E2] 5i — > 62) are also forbidden transitions with Av = in the 7-unstable 
limit, they are not very small in Fig. 20 and indicate a significant deviation from the 
7-unstable limit. 

The spectroscopic quadrupole moments Q calculated in the (1+3)D and the 
(2+3)D models are compared in Fig. 21. It is seen again that the two calculations 
show the same qualitative features: both calculations yield the positive (negative) 
sign for the spectroscopic quadrupole moments of the yrast (yrare) states to indicate 
an oblatelike (prolatelike) character. Quantitatively, in the (2+3)D model, the yrast 
Q value increases with angular momentum more significantly and the absolute values 
of the yrare Q moments are slightly smaller than those of the (1+3)D model. As we 
discussed above, in both calculations, the Q value of the yrare states approaches zero 
with the angular momentum increasing due to the cancellation mechanism associated 
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Fig. 21. The same as Fig. 20 but for the spectroscopic quadrupole moments in units of the intrinsic 
quadrupole moment, eQo = 3e/\/57r^-Ro/3o. The results of calculation in the (1+3)D model and 
in the (2+3)D model are displayed as functions of angular momentum / in the left and right 
panels, respectively. The values for the yrast (yrare) states are shown by filled (open) squares. 



with the growth of the two peak structure in the collective wave functions. In spite of 
such deviations from a simple picture, we can see in Fig. 21 some qualitative features 
characteristic to the oblate-prolate shape coexistence. 

We conclude that the excitation spectrum and the properties of the quadrupole 
transitions and moments exhibited in Figs. 20 and 21 can be regarded as those 
characteristic to an intermediate situation between the well-developed oblate-prolate 
shape coexistence and the 7-unstable limit. 

§6. Concluding Remarks 

From a viewpoint of the oblate-prolate symmetry and its breaking, we have pro- 
posed a simple (1+3)D model capable of describing the coupled motion of the large- 
amplitude shape fluctuation in the 7-degree of freedom and the three-dimensional 
rotation. Using this model, we have made a systematic investigation of the oblate- 
prolate shape coexistence phenomena and their relationships to other classes of low- 
frequency quadrupole modes of excitation, including particular cases described by 
the 7-unstable model and the rigid triaxial rotor model. We have also adopted the 
(2+3)D model to check the validity of freezing the (3 degree of freedom in the (1+3)D 
model. We have obtained a number of interesting suggestions for the properties of 
low-lying states that are characteristic to an intermediate situation between the well- 
developed oblate-prolate shape-coexistence and 7-unstable limits. In particular, 1) 
the relative energies of the excited + states can be indicators of the barrier height 
of the collective potential. 2) Specific E2 transition probabilities are sensitive to the 
oblate-prolate symmetry breaking. 3) Nuclear rotation can assist the localization of 
the collective wave functions in the (/3,7) deformation space. However, even if the 
rotation-assisted localization is realized in the yrast band, it is not necessarily in the 
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yrare band: the two-peak structure may develop in the yrare band. 
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